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During its first year of data taking, the Large Area Telescope (LAT) onboard the Fermi Gamma- 
Ray Space Telescope has collected a large sample of high-energy cosmic-ray electrons and positrons 
(CREs). We present the results of a directional analysis of the CRE events, in which we searched 
for a flux excess correlated with the direction of the Sun. Two different and complementary analysis 
approaches were implemented, and neither yielded evidence of a significant CRE flux excess from 
the Sun. We derive upper limits on the CRE flux from the Sun's direction, and use these bounds 
to constrain two classes of dark matter models which predict a solar CRE flux: (1) models in which 
dark matter annihilates to CREs via a light intermediate state, and (2) inelastic dark matter models 
in which dark matter annihilates to CREs. 
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I. INTRODUCTION 

In the last decades the searches for a dark matter 
(DM) signal from the Sun were performed looking for 
possible excesses of neutrinos or gamma-rays associated 
with the Sun's direction. However, as it was noted in 
Ref. pj], several DM models that have been recently 
developed to explain various experimental results also 
imply an associated solar flux of high-energy cosmic- 
ray electrons and positrons (CREs). On the other 
hand, no known astrophysical mechanisms are expected 
to generate a significant high-energy CRE (> 100 GeV) 
excess associated with the Sun. 

A class of models in which DM annihilates to CREs 
through a new light intermediate state </» @, Q has been 
considered to explain the excesses in local CRE fluxes 
reported by PAMELA Q, ATIC Q, and Fermi I, Q. 
In these scenarios DM particles captured by the Sun 
through elastic scattering interactions would annihilate 
to 4> pairs in the Sun's core, and if the (j) could escape 
the surface of the Sun before decaying to CREs, these 
models could produce an observable CRE flux. 

Another class of models in which DM scatters off 
of nucleons predominantly via inelastic scattering has 
been proposed as a means of reconciling the results of 
DAMA and DAMA/LIBRA @Ji with CDMS-II [lo|,[ll| 
and other experiments (e.g., [1^, [Ij]; see also [ij] for 
a comprehensive discussion of experimental constraints). 
If DM is captured by the Sun only through inelastic 
scattering (IDM), this could lead to a non- negligible 
fraction of DM annihilating outside of the Sun's surface. 
For models in which iDM annihilates to CREs, an 
observable flux at energies above a few tens of GeV could 
be produced. 
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During its first year of operation, the Lar ge Area 
Telescope (LAT) onboard the Fermi satellite [l^ has 
collected a substantial number of CRE events, which has 
allowed a precise measurement of the energy spectrum 
over a broad energy range from a few GeV up to 1 TeV p, 
'7|. Furthermore, a directional analysis of the high-energy 
CRE events was performed in the Galactic reference 
frame [16], and showed no evidence of anisotropics . 

In this paper we use the high-energy CRE data set 
to search for flux variations correlated with the Sun's 
direction. Since the Sun is moving with respect to the 
Galactic reference frame, the previously-reported absence 
of anisotropics in the CRE flux observed in the Galactic 
frame does not necessarily imply a negative result. 



II. DATA SELECTION 

The Fermi LAT is a pair-conversion telescope designed 
to detect gamma rays in the energy range from 20 MeV to 
more than 300 GeV. A full description of the apparatus 
is given in 15]. Even though it is a photon detector, it 
has been demonstrated that the LAT is also an excellent 
CRE detector [Q, [7i, JJi]. For this analysis we used the 
CRE data sample collected by the LAT during its first 
year of operation, starting from August 4, 2008. The 
event _selection was performed in the same way as in 
Ref. 16]; approximately 1.35 x 10^ CRE events with 
energies larger than 60 GeV passed the selection cuts. As 
discussed in Ref. [1^ , the energy threshold of 60 GeV was 
chosen because it is higher than the geomagnetic cutoff 
in any part of Fermi's orbit. 

Unlike gamma rays, CREs are deflected by interactions 
with magnetic fields encountered during their propaga- 
tion in interstellar space. In particular, CREs coming 
from the Sun are deflected by both the Sun's and the 
Earth's magnetic fields. 

Geomagnetic effects on CREs have been studied using 
a code that reconstructs the trajectories of charged par- 
ticles in the Earth's magnetic field based on the Interna- 



tional Geomagnetic Reference Field (IGRF) model [17| . 
Since the Fermi LAT cannot measure the sign of the elec- 
tric charge, we associated both an electron and a positron 
track with each CRE event detected by the LAT. Each 
track starts from the detection point with the same en- 
ergy of the event and with a direction opposite to that of 
the event, and ends at a very large distance (larger than 
100 Earth-radii from the Earth's center). 

The distribution of deflection angles at different ener- 
gies was analyzed. The simulation demonstrated that, 
at energies above 20 GeV, 90% of the particles are de- 
flected with respect to the original direction within an 
angle ^90% given by the approximate formula: 



2.8° 



E{ TeV) 



(1) 



where E is the particle energy. Hence, due to the 
geomagnetic field, the reconstructed directions of CREs 
with energies above 100 GeV detected by the LAT and 
coming from any given direction of the sky will be spread 
over a cone with an angular radius of about 30° centered 
on the original incoming direction. 

The directions of incoming CREs are also affected by 
the Heliospheric Magnetic Field. A detailed study of 
its effects on CREs is beyond the scope of this work, 
however in Ref. 18], it was shown that CREs with 
energies of several hundreds of GeV can travel through 
the center of the solar system without experiencing 
significant defiections. 

CREs travelling in the Solar System may also suffer 
energy losses, mainly due to the Inverse Compton (IC) 
scatterings on the photons emitted by the Sun and to the 
synchrotron radiation (SR) emitted in the interactions 
with the Heliospheric Magnetic Field. To study the 
energy loss processes of CREs travelling from the Sun 
to the Earth we implemented a simple toy model, in 
which we assumed that CREs propagate from the Sun's 
surface to the Earth in straight lines and with velocity 
c. Following Ref. [T^, we assumed that the Sun can be 
modeled as a black body with a temperature of 5777 K 
and with a photon density given by: 



density evaluated from eq. [2l (Jt is the Thomson cross 
section and jk is given by: 



Ik 



iy/birLeC^ 



(4) 



where rUe is the electron mass, ks is the Boltzmann's 
constant and T = 5777K is the temperature of the 
Sun's surface. The evaluation of the SR energy loss rate 
is not easy, because the structure of the Heliospheric 
Magnetic Field is rather complex [2l|. However, as a 
first approximation, we assumed that the strength of the 
Heliospheric Magnetic Field drops from the Sun's surface 
as: 
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where Bq = 1 gauss is the strength of the field on the 
Sun's surface. We did not include the contribution of the 
Geomagnetic field to synchrotron energy losses because, 
even though the field strength at the LAT altitude is 
of the order of 1 gauss, the path length of CREs in the 
Geomagnetic field is of the order of a few Earth radii, 
which is negligible with respect to the path length in the 
Heliospheric Magnetic field, which is of the order of a few 
solar radii. 

The synchrotron energy loss rate of CREs was then 
calculated as [221: 
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where Wb is the magnetic field energy density that 
includes only the contribution from the Heliospheric 
Magnetic Field and, in our model, results to be negliglible 
with respect to the IC energy loss rate. Using eqs. [3] 
and [HI we calculated that CREs in the energy range from 
60 GeV to 1 TeV travelling from the Sun to the Earth lose 
no more than 2% of their initial energy. Therefore, in the 
calculations of the following sections, we will neglect all 
energy loss processes. 



Nphie,r) = 0.5n;,b(e) 
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where nbb{^) is the blackbody photon energy density 
(Planck's equation), i?© is the solar radius and r is the 
distance from the center of the Sun. The IC energy loss 
rate of CREs was then evaluated as in Ref. 1201 as: 



III. DATA ANALYSIS AND RESULTS 

To study the CRE flux from the Sun's direction and 
to search for variations with respect to the average flux, 
we implemented two complementary analysis approaches: 
(i) flux asymmetry analysis and (ii) comparison of the 
solar flux with the isotropic flux. 
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where /?c and 7 are respectively the velocity and the 
Lorentz factor of the CRE, W is the photon energy 



A. Flux asymmetry studies 

This approach compares the CRE flux from the Sun 
with the flux from a fake source (fake Sun) placed in the 
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FIG. 1: Differential flux asymmetry between real and fake Sun evaluated in cones with angular radii AG = 30° (top left panel), 
45° (top right panel), 60° (bottom left panel) and 90° (bottom right panel). The fluxes are multiplied by E^ (the energy values 
correspond to the bin centers) since the energy spectrum of CREs is approximately proportional to E^'' in this energy range. 
Only statistical error bars are shown. 



Angular radius 


Maximum deviation (amax) 


P{\a„,ax\) 


P{W\ > \(T^ax\) 


30° 


2.690 


0.007 


0.113 


45° 


-2.542 


0.011 


0.171 


60° 


-2.806 


0.005 


0.082 


90° 


-2.947 


0.003 


0.050 



TABLE I: For each cone used for the flux asymmetry analysis the maximum deviations (either positive or negative) from the 
null value are shown and the corresponding probabilities of observing larger values in the hypothesis of null flux asymmetry. 
The last column shows the probability of finding at least one energy bin with a larger flux asymmetry than the maximum 
observed value. 



sky position opposite to that of the Sun. To perform our 
analyses, we chose a custom reference frame derived from 
ecUptic coordinates. The ecliptic coordinates associated 
with each CRE event were evaluated from equatorial 
coordinates using the formulae in Ref. 23]. Indicating 
with (A, /3) the pair of ecliptic coordinates (longitude 
and latitude, respectively) associated with any given 
direction, the Sun's direction will always lie in the plane 
/3 = 0. In fact, since the Sun is moving eastwards along 
the path of the ecliptic, its ecliptic latitude will always 



be zero by definition, while its ecliptic longitude will 
always increase, describing a complete 360° cycle in one 
year [23J . In our custom reference frame, the coordinates 
associated with each direction are defined as: 



A' = A — Xsu 

13' ^P 



(7) 



where Xsun is the echptic longitude of the Sun, evaluated 
from the Sun's ephemeris using a software interfaced to 
the JPL libraries [2J]. In this reference frame, the Sun's 
coordinates wiU always be {X'g^,^ — 0°,/35„„ = 0°). On 
the other hand, the fake Sun will always be located at 
the coordinates (A},,, sun = 180°, /3},,, sun = 0°). 

Due to the geomagnetic field's effects on CRE trajecto- 
ries described in ^TTJ we consider the fluxes from extended 
sky regions centered on the Sun (and on the fake Sun). 
In particular, we compare the CRE fluxes from directions 
within cones of angular radii A0, centered on the posi- 
tion of the Sun and the fake Sun. According to Eq.[Tl 90% 
of CREs with energies of 100 GeV are deflected within a 
cone of about 30° angular radius, and so we chose this 
value as the minimum angular radius of the sky regions 
to be investigated because the DM models discussed in 
Ref. [l| predict a CRE fiux excess from the Sun in the 
energy range above 100 GeV. 

To measure the fluxes from different sky regions, we 
flrst divided the sky into a grid of pixels, then evaluated 
the CRE fluxes from individual pixels (each pixel was 
treated as a point source), and flnally integrated the 
fluxes from the pixels belonging to the selected sky 
regions. We used the HEALPix [25] pixelization scheme, 
and divided the sky into 12288 equal-area pixels, each 
covering a solid angle of about 10~^ sr. The CRE 
differential fluxes from individual pixels were evaluated 
according to the following equation: 
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where d^i{E)/dE is the differential CRE flux (expressed 
in particles per unit energy, unit area and unit time) 
in the energy interval [E,E + AE] from the ith pixel, 
Ni{E) is the number of observed CRE events from the 
ith pixel with energies between E and E -f AE, c{E) 
is the residual contamination (the contamination values 
are reported in Ref. Q) and £i{E) is the exposure of 
the ith pixel, which is calculated taking into account the 
effective area of the instrument, and the live time of the 
ith pixel. The dependence of the effective area on the 
CRE direction in the instrument, expressed in terms of 
the off-axis and azimuth angles 9 and 0, is also taken into 
account in the calculation. 

The CRE flux from a cone of angular radius AG 
centered on the Sun is then given by: 



d^Sun{E\Ae) 
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where ROI{AQ) denotes the set of pixels (region of 
interest) at an angular distance less than A 9 from the 
Sun. The flux from the fake Sun is evaluated in a similar 
way. The flux asymmetry can then be evaluated as: 



dA^iE\AQ) d'PsuniElAQ) d<^Fake SuniE\AQ) 
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The variable dA^{E\Ae)/dE defined in Eq. [TU] is the 
difference between the CRE flux from the Sun and the 
fake Sun; the flux of the fake Sun is assumed to be 
representative of the average CRE flux across the sky. 
Positive (negative) values oi dA,s>{E\AQ) /dE indicate an 
excess (deficit) of CREs from the Sun. We emphasize 
that this approach relies on the assumption that the flux 
from the fake Sun region is representative of the average 
CRE flux. 

In Fig. [T] the differential CRE flux asymmetries 
dA^{E\AQ)/dE between the real and the fake Sun are 
shown for four different ROIs, with angular radii of 30°, 
45°, 60° and 90°. No signiflcant CRE flux excesses or 
deficits from the Sun are observed at any energy. In 
the plots of Fig. [1] only statistical error bars are shown. 
As pointed out in Ref. [7|, the main source of system- 
atic uncertainties in the evaluation of CRE fluxes is the 
imperfect knowledge of the detector's effective area. As- 
suming that it is affected only by a normalization er- 
ror, when calculating the error on the flux differences 
dA^{E\AQ)/dE, the contribution from the normaliza- 
tion error will be proportional to \dA^{E\AQ) / dE\ (see 
the discussion in Ref. [26]), and therefore it will be neg- 
ligible with respect to the statistical error. 

Assuming that the measured flux asymmetries in each 
energy bin behave as Gaussian random variables, we 
expressed the excesses and deflcits (with respect to the 
hypothesis of a null flux asymmetry) in units of a {a is the 
statistical error associated with each measurement), and 
evaluated the corresponding probabilities of measuring 
larger excesses or deficits assuming the null hypothesis. 
Table |T] shows, for each value of the angular radius 
AQ, the maximum observed deviations from the null 
flux asymmetry in units of cr, and the corresponding 
probabilities of measuring larger flux asymmetries in 
the null hypothesis. As shown in Table U the flux 
asymmetries in all of the ROIs are always within Sa of 
zero. The last column of Table U shows the probabilities 
of flnding, in each ROI, at least one energy bin with a 
flux asymmetry larger than the maximum observed value. 
The probabilities were calculated assuming that the flux 
asymmetries measured in each of the 17 energy bins used 
in our analysis are uncorrelated. The calculations were 
performed taking only statistical errors into account; 
if systematic errors were also taken into account, the 
signiflcance of the deviations of the flux asymmetries 
from zero would be smaller. 

The same analysis was repeated for integral fluxes 
above various energy thresholds and again no evidence 
of flux asymmetries was found. 



1. Evaluation of statistical upper limits on the C RE flux 
asymmetry 



The count difference is finally converted into a flux 
asymmetry according to the following equation: 



The previous analysis did not provide any evidence 
of a CRE flux excess from the Sun with respect to 
the fake Sun, so we set statistical upper limits on this 
signal by following the approach outlined in Ref. f27'| 
(pp. 136-139). In each energy bin, the measured flux 
asymmetry dA^{E\lS.<d) / dE can be seen as a realization 
of a Gaussian random variable. Assuming the hypothesis 
of a CRE flux excess from the Sun, its expectation value 
must be non-negative: 



' dE 



>0 



(11) 



To set upper limits on {dA^{E\l:^Q) / dE) we imple- 
mented the Bayesian method described in Ref. [23, as- 
suming a uniform prior density. The a of the Gaussian 
probability distribution function associated with each 
measurement of dA^{E\/\Q) / dE is the statistical error 
associated with the measurement. 

Fig. [5] shows the statistical upper limits on CRE flux 
asymmetries at different confidence levels, considering 
different cones centered on the Sun with angular radii 
ranging from 30° to 90°. In Table HIl the values of the 
upper limits at 95% confidence level for the flux asym- 
metry (iyl$(£'|A8)/(i£' in the different ROIs are summa- 
rized. The calculated upper limits are also expressed in 
terms of fractions of the CRE flux from the region of the 
fake Sun. 



B. Comparison with an isotropic flux 

The second approach used in this analysis is based on 
the event-shuffling technique employed in Ref. [16], which 
was used to build a simulated sample of isotropic CREs 
starting from the real events. Simulated events are built 
by randomly coupling the arrival times and the arrival 
directions (in local instrument coordinates) of real events. 
The simulated eventsample used for this analysis is the 
same used in Ref. jl_6| and is 100 times larger than the 
real one. 

In this case, given the angular radius AG of a cone 
centered on the Sun, we evaluated the count differences 
between real and simulated CREs as: 



dA'^{E\^Q) 
dE 



1 Aiv(£;|Ae)(i 
~Ke 



c{E)) 



£{E\/^Q) 



(13) 



where c{E) is the residual contamination and £{E\^Q) is 
the exposure of the selected sky region, which is evaluated 
from the instrument effective area and the live times of 
the pixels belonging to that sky region. 

We emphasize that the two approaches used in this 
work are complementary, but not fully equivalent. In 
the first case the variable dA,^{E\/S.<d) / dE is built using 
real events from different directions (real Sun and fake 
Sun). On the other hand, in the second case, the variable 
dA'^{E\/S.Q) / dE is built using real and simulated events 
from the same region of the sky. In both cases the goal 
of the analysis is to compare the CRE flux from the Sun 
with the average CRE flux. In the first case, the reference 
flux is evaluated looking at real events from the fake Sun 
region, while in the second case it is evaluated simulating 
an isotropic flux and looking at simulated events from the 
Sun region. 

The second approach, however, excludes potential 
systematic biases when calculating flux differences. In 
particular, to evaluate the flux from a given region of 
the sky requires knowledge of the exposure, which in 
turn depends on the effective area of the detector and 
on the observation live time. The effective area is 
calculated from Monte Carlo simulations and thus could 
be affected by systematics such as variations correlated 
with time and spacecraft position or miscalculations of its 
dependence on instrument coordinates. When evaluating 
the flux asymmetry dA,^{E\/S.<d) / dE according to Eg.lTUl 
the systematic uncertainties involved in the evaluation of 
the two terms could be different, and the result could be 
biased. On the other hand, when evaluating the flux 
difference dA'i^{E\IS.<d) / dE from Eq. [T31 inaccuracies in 
the effective area calculation can only result in a scale 
error on the flux difference. 

To determine whether the real counts differ signif- 
icantly from the simulated ones, we performed a hy- 
pothesis test following the prescriptions of Ref. j28| . 
Denoting with Nreai{E\I^<d) and Nsim{E\/^Q) the real 
and the simulated counts in the energy bin _E, the 
null hypothesis for this analysis is that Nreai{E\/S.<d) — 
a{E)Nst,niE\Ae). Following Ref. [H], we evaluated the 
significance in each energy bin as: 



AiV(^|Ae) = Nreai{E\Ae) - a{E)Nsun{E\AQ) (12) 

where Nreai{E\AQ) and Nsim{E\AQ) are respectively 
the number of CRE events in the real and simulated 
data sets with energy E and arrival directions within 
the selected cone. The parameter a{E) in Eq. [12] is a 
normalization factor. If Nreai{E) and Nsim{E) are the 
total numbers of real and simulated events with energy 
E, then aiE) = NreaiiE)/Ns^miE) M- 



S iE\AQ) ^ ±V2 I^Nreai{E\ AQ)\n 



Nreal{E\A&) 



l + a{E) 



Nreal{E\Ae) + N,,„,{E\Ae) 

N,,rr.iE\Ae) 



a{E) 
Ns„n{E\AQ) 



In 



{l + a{E)) 



Nreal{E\AQ) + N,,^{E\AQ)_ 



1/2 



(14) 
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FIG. 2: Statistical upper limits at confidence levels of 68% (dotted lines), 95% (dashed lines) and 99% (continuous lines) for 
the CRE flux asymmetry between real and fake Sun, evaluated in cones with angular radii A0 = 30° (top left panel), 45° (top 
right panel), 60° (bottom left panel) and 90° (bottom right panel). 



with the convention of choosing the + sign if 
Nreai{E\Ae) > a{E)Ns^m{E\Ae) and the - sign if 
Nreai{E\AQ) < a{E)Ns^„^iE\Ae) . 

The significance values evaluated from Eq. [14] can be 
converted into probability values. In particular, since 
S^ is a random variable following a x^ distribution with 
1 degree of freedom [2a|, one can easily evaluate the 
probability of observing a value of S^ larger than the one 
observed. In Table Hill the maximum deviations from the 
null flux asymmetries are shown for each ROI used for our 
analysis, together with the corresponding probabilities of 
finding larger deviations. The last column of the table 
shows the probabilities of finding, in each ROI, at least 
one energy bin with a flux asymmetry larger than the 
maximum observed value. Again, these probabilities 
were calculated assuming that the flux asymmetries 
measured in each of the 17 energy bins used for our 
analysis are uncorrelated. The observed deviations from 
the null flux asymmetries are statistically insignificant. 

The count differences were converted into flux differ- 
ences according to Eq. 1131 In Fig. |3] the asymmetry 
variable dA'^{E\AQ) / dE between the real and simulated 



fluxes from the Sun is shown for the four cones with an- 
gular radii of 30°, 45°, 60° and 90°. Again, no evidence 
of a CRE signal from the Sun is observed. Similar results 
are obtained when integral fluxes are analyzed. 



Evaluatio 



<n of statistical upper limits on the CRE flux 



from the Sun 



As in ijIII Al our analysis in this case does not find evi- 
dence of a CRE signal from the Sun, and so we set statis- 
tical upper limits on solar CRE fluxes. Given a cone with 
angular radius AG centered on the Sun, the observed 
counts Nreai{E\AQ) in the energy bin E can be seen as 
a realization of a Poisson random variable. Assuming the 
hypothesis of a CRE signal from the Sun, these counts 
will be the sum of a signal contribution, Nsun{E\AQ), 
plus a background contribution, Nhkg{E\AQ): 



NrealiE\AQ) = NsuniE\AQ) + NbkgiE\Ae). (15) 
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TABLE II: Statistical upper limits at 95% confidence level on the CRE flux asymmetries between the real and the fake Sun, 
evaluated in cones of angular radii Ae — 30° , 45° , 60° , 90° . The upper limits are also expressed in terms of fractions of the 
CRE flux from the fake Sun. 



Angular radius 


Maximum deviation (Smax) 


Pi\Smax\) 


P{\S\ > \Smax\) 


30° 


1.925 


0.054 


0.611 


45° 


0.419 


0.675 


1.000 


60° 


-0.654 


0.513 


1.000 


90° 


1.026 


0.305 


0.998 



TABLE III: For each cone used for the flux asymmetry analysis the maximum deviation (either positive or negative) in terms 
of significance from the null value are shown with the corresponding probability of finding a larger significance value. The last 
column shows the probability of finding at least one energy bin a with larger flux asymmetry than the maximum observed 
value. 



Both Nsun{E\AQ) and Niikg{E\AQ) can be seen as Pois- 
son random variables. The only information available 
about iV5„„(£'|A8) is that, in the hypothesis of a CRE 
signal from the Sun, its average value (which is unknown) 
must be non-negative: 



{Nsun{E\AQ)) > 0. 



(16) 



On the other hand, the average value of Nbkg{E\AQ) can 
be estimated from the randomized data sets, and is given 
by: 



{NbkgiE\AQ)) ^ aiE)Ns^ 



{E\A&) (17) 

total real and 



where a{E) is the ratio between the 
simulated events in the energy bin E. 

The goal of our analysis is to evaluate an upper limit on 
the average value (A^5„„(£'|A8)), that will be converted 



into an upper limit on the CRE flux from the Sun after 
properly taking into account the detector acceptance and 
the live time. For our calculation we implemented the 
Bayesian method with the assumption of a uniform prior 
density. The mathematical details of the method can be 
found in Ref. [13 (pp. 139-142). 



Fig. m shows the statistical upper limits at different 
confidence levels on the CRE fluxes from different cones 
centered on the Sun and with angular radii ranging from 
30° to 90°. The results are consistent with those shown in 
Fig.m In Table HV] the values of the upper limits at 95% 
confidence level for the flux asymmetry dA'^{E\AQ)/dE 
in the different ROIs are summarized. The calculated 
upper limits are also expressed in terms of fractional 
excess with respect to the expected isotropic flux. 
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FIG. 3: Differential fiux asymmetry between real and simulated events from the Sun evaluated in cones with angular radii 
AO — 30° (top left panel), 45° (top right panel), 60° (bottom left panel) and 90° (bottom right panel). Only statistical error 
bars are shown. 



2. Spherical harmonics analysis 

The previous analyses excluded variations of the CRE 
flux correlated from the Sun's direction. However, due 
to the effects of the heliospheric magnetic field and the 
geomagnetic field, a CRE signal from the Sun could 
produce a flux excess from a direction shifted with respect 
to the Sun's position. Also, since CREs from the Sun are 
expected to be spread over a cone with a finite angular 
radius, an excess of CREs from the Sun could induce an 
anisotropy on a large angular scale. 

To investigate this possibility we implemented a more 
general analysis method, based on spherical harmonics 
analysis of a fluctuation sky map. This method was 
applied in Ref. [1^] to search for anisotropics in the CRE 
flux in the Galactic reference frame, while in this paper 
we adopted the custom coordinates in Eq. [7] derived from 
the ecliptic reference frame. 

The fluctuation sky map was built starting from the 
real sky map and from the simulated one, which was 
generated using the randomized data sets. The analysis 
was performed on sky maps of the counts integrated 



above a given energy in order to retain sufficient statistics 
in each energy bin. The fluctuation in the ith pixel is: 



/.(> E) = 



Ni^real{> E) - aNi^,i„^{> E) 



aNi_ 



.{>E) 



(18) 



The fluctuation sky map is then expanded in the basis 
of spherical harmonics to obtain the set of coefficients 
aim- The coefficients of the angular power spectrum are 
given by the variance of the 2/ + 1 aim coefficients at each 
multipole: 



Ci = 



21 + 1 ^^ 



I 

\~^ I |2 

aim] 



(19) 



m— - 



Each coefficient Ci characterizes the intensity of the 
fluctuations on an angular scale of ~ 180°//. In the case 
of an isotropic flux, each of the coefficients C; can be 
seen as a random variable with a true value (white noise) 
given by: 
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FIG. 4: Statistical upper limits at confidence levels of of 68% (dotted lines), 95% (dashed lines) and 99% (continuous lines) 
for the ORE fluxes from the Sun, evaluated in cones with angular radii AO — 30° (top left panel), 45° (top right panel), 60° 
(bottom left panel) and 90° (bottom right panel). 



IV. SOLAR CRE FLUXES FROM DARK 
MATTER 



Ci^C 



N 



An 



(20) 



where N is the total number of observed events. The 
confidence intervals for the Ci can be evaluated by noting 
that the random variable {21 + l)Ci/Ci follows a X2i+i 
distribution. 

In Fig. [S]the angular power spectra after the subtrac- 
tion of the white noise contribution Cjv are shown for four 
different minimum energies (60 GeV, 100 GeV, 200 GeV 
and 500 GeV). The curves show the 3ct and 5cr proba- 
bility intervals assuming the hypothesis of an isotropic 
CRE flux. All the data points lie within the 3a interval, 
indicating that the measurements are consistent with the 
hypothesis of an isotropic CRE flux. Hence, we conclude 
that no preferred CRE arrival directions are observed. 



We now determine constraints on DM model parame- 
ters by comparing our upper limits on solar CRE fluxes 
to the predicted fluxes of the two DM annihilation sce- 
narios considered in Ref. [ij: (1) capture of DM particles 
by the Sun via elastic scattering interactions and subse- 
quent annihilation to e through an intermediate state 
(/), and (2) capture of DM particles by the Sun via inelas- 
tic scattering interactions and subsequent annihilation of 
the captured DM particles outside the Sun directly to 



Dark matter annihilation through an 
intermediate state 



In this case we assume the standard scenario for WIMP 
capture by the Sun, namely that DM particles x are cap- 
tured by the Sun through elastic scattering interactions 
and then continue to lose energy through subsequent 
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TABLE IV: Statistical upper limits at 95% confidence level on the CRE flux asymmetries between the real and the simulated 
Sun, generated by the event shuffling technique, evaluated in cones of angular radii Ae — 30°, 45°, 60°, 90°. The upper limits 
are also expressed in terms of fractions of the isotropic CRE flux from the simulated Sun. 



scatterings, eventually thermalizing and sinking to the 
core where they annihilate. In general, the only annihi- 
lation products which can escape the Sun are neutrinos; 
photons and charged particle final states are trapped by 
interactions with the dense matter in the Sun. However, 
recently scenarios have been proposed in which DM par- 
ticles annihilate into a light intermediate state 0, i.e., 
XX ~^ 4'4'i with the (p subsequently decaying to standard 
model particles; these models have been suggested to pro- 
vide a means of explaining an excess in the CRE spec- 
trum reported by ATIC and Fermi and in the positron 
fraction reported by PAMELA by DM annihilation or 
decay |29l - l35j . For the case considered in Ref. pj], the 
are assumed to be able to escape the Sun without fur- 
ther interactions, with each decaying to an e pair. If 
this decay happens outside the surface of the Sun, the e^ 
could reach the Earth and may be detectable in the form 
of an observed excess of CREs from the direction of the 
Sun. 

The DM particles are assumed to annihilate at rest in 
the core of the Sun, so in the lab frame the energy of the 
0, E^, is equal to the DM particle mass m^. We assume 
(/) to be a light scalar such that rn^ <C m^, hence the 
are relativistic. The energy of the (j) is described by the 
parameter /3ci = Vd/c where Vc\ is the relative velocity of 
the lab frame and the </> rest frame (hereafter, the CM 
frame). The </> are assumed to lose a negligible amount 
of energy exiting the Sun. A decays into an e pair 
with an isotropic angular distribution in the CM frame. 
Both the e+ and e" have the same energy in the CM 



frame, parameterized by /3jc = I'jc/c = i/l — {Anil/rrii), 

where me is the electron mass and Wjc is the velocity of 
particle j (the electron or positron) in the CM frame. 
However, in the lab frame the e^ are boosted and so the 
angular distribution is no longer isotropic, and the energy 
of an e in the lab frame depends on the angle at which 
it is emitted relative to the direction in which the 4> is 
traveling. In the lab frame the angle at which the e is 
emitted is denoted ^lab; in the CM frame it is 9 cm- The 
e are assumed to travel in straight lines and to suffer no 
energy losses before reaching the detector, i.e., the effects 
of magnetic fields are assumed to be negligible. 

The flux of e per time, area, and energy per cosine 
of the detector angle from the direction 0det is given by 
integrating the differential rate of decay of (j) along the 
line of sight in the direction of i 



dN 



dt dA d cos 0dctc?£'dct 
°°^j^ dN dT 

dVdt dcos^dct 



7dct, 

^dctj-Edot) = 

<5(^dct - E), 



(21) 



where R is the distance from the detector in the line-of- 
sight direction defined by ^dct- At each 0dct we exclude 
the contribution from (p decays occurring in the range of 
R within the surface of the Sun. 

The first term in Eq. [21] is the rate of production of 
e^ (equal to twice the rate of 4> decays) per volume at a 
volume element a distance r from the center of the Sun, 
and is given by 
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FIG. 5: Angular power spectra for different minimum energies: 60 GeV (top left panel), 100 GeV (top right panel), 200 GeV 
(bottom left panel), 500 GeV (bottom right panel). The points show the quantities Ci — Cn- The dashed lines and the 
continuous lines show respectively the 3a and 5a intervals for the probability distribution of the white noise. 



dN 
dVdt 



(r(0dct,i?)) 






(22) 



where C© is the capture rate of DM particles in the Sun 
and the characteristic decay length L = 7cict, where 
7ci = 7(/3ci) with 7(/3) — l/-\/l — /3'^ and r is the lifetime 
of the (j). Equilibrium is assumed, i.e., for every two 
DM particles that are captured, two annihilate, so C© 
also represents the rate of production of (p particles by 
annihilation at the center of the Sun. For the range of 
scattering cross-sections considered here, the capture rate 
is in general sufficiently large that equilibrium is a valid 
assumption, although we discuss this issue in greater 
detail when presenting our constraints. 

We consider separately the cases of solar capture by 
spin-independent scattering and spin-dependent scatter- 
ing. Using DarkSUSY [36| , we calculate the capture rate 
via elastic scattering C© as a function of m^ . A Maxwell- 
Boltzmann velocity distribution is assumed, with the so- 
lar velocity relative to the DM rest frame u© = 220 km/s, 
the DM velocity dispersion v = 270 km/s, and the lo- 



cal DM density pDM — CSGeV/cm"^. For the case of 
spin-independent scattering we set the spin-dependent 
scattering cross-section ctsd = Ocm^ and calculate the 
capture rates for ctsi — lO^'^'^cm^; for spin-dependent 
scattering we set csi — Ocm^ and calculate the capture 
rates for ctsd = 10"'*" cm^. The reference values of crsi 
and (TsD roughly correspond to the current experimental 
upper limits on these parameters. For the range of pa- 
rameters considered here the capture rate scales linearly 
with cTgD and ctsi. 

The second term in Eq. [5T] is the angular distribution 
of the e^ from the decay of a 0, as observed in the lab 
frame, and expressed in terms of detector angle. 



dV 



dV 



d cos 9, 



dct 



d cos 9r 



d cos 6c 



d cos 9 



lab 



d COS 9 



lab 



d COS 9, 



dct 



(23) 



In the rest frame of the </> the decays are isotropic, so, 
after integrating over the azimuthal angle we can write 



dV 



d cos 9c 



-1/2. 



(24) 
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The transformation between the CM angle and lab angle 
is given by [S^l 



d cos 9r 



d cos 6'lab 



[jlja + cos e^^f + sin^ 9cru?^^ 
|7ci(l + acos6'cm)| 



, (25) 



with 7ci — 7(/3ci) and a — /3ci//?jc- The lab and detector 
angles are related by 



71ab = Cdot 



sni 



J /_R sin 6*^01 



(26) 



which gives 



d cos 6'lab _ (I-D0 - J?cos(6ldct)| + ficos(6'dct))^ , . 



d cos 6'dGt 



Dq -i?cos(6'dct) 



The delta function in Eq. [5T] enforces that the energy 
observed at the detector is equal to the energy of the 
emitted e^ boosted to the lab frame, 
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FIG. 6: Constraints on DM annihilation to e^e~ via an 
intermediate state, from solar CRE flux upper limits. Solar 
capture of DM is assumed to take place via spin-independent 
scattering. The constraints obtained for three values of the 
decay length L of the intermediate state are shown. Models 
above the curves exceed the solar CRE flux upper limit at 
95% CL for a 30° ROI centered on the Sun. 



-E'(^cin) = 2 701^0(1 + ^cl COS( 



0- 



(28) 



Note that because the energy in the lab frame depends 
only on 0cm, and because ^lab is determined by 0cm, 
fixing i?dct corresponds to selecting only CREs emitted 
at the corresponding 6'lab • For a specified 6'dct, the 6'lab of 
particles observed along the line-of-sight R varies, hence 
the observed energy of CREs emitted from a point along 
the line-of-sight is a function of R, i.e., E^ctiR)- We 
rewrite the delta function in Eq. [5T] as the composition 



5{Ei,t - E{R)) 



6iR - Rq) 

di^(^") 



(29) 



and then perform the integration over R. The parameter 
Rq is the value of R along the line-of-sight in the direction 
0dot where 0iab takes the value required to generate CREs 
with a given ii^dct- 

We evaluate the CRE flux within a ROI of 30° centered 
on the Sun, and fix the value of m^ — 1 GeV. We 
calculate limits for three values of the decay length 
L = 5AU, lAU, and 0.1 AU. Decreasing L increases 
the observed CRE flux by condensing the region within 
which most (j) decay. However, we emphasize that even 
for as large a decay length as L = 5 AU, the signal 
in the energy range used in this analysis is strongly 
peaked in the direction of the Sun and extends only a 
few degrees at most. Since the (j) in this scenario are 
relativistic, in the lab frame the emitted e^ are boosted 
along the direction the </> is moving, and so only (j) exiting 
the Sun very close to the direction of the detector will 
produce decay products with large enough 0iab to reach 
the detector. In particular, for the e to have sufficient 
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FIG. 7: Constraints on DM parameters for annihilation to 
e^e~ via an intermediate state as in Fig. [6l except assuming 
solar capture by spin-dependent scattering. 



energy to fall within the energy range of this analysis, 
a significant fraction of the energy must be deposited 
into the e^ that reach the detector. This only occurs 
for e^ emitted with very small 0iab. This also leads 
to an energy dependence of the angular signal: for a 
given DM scenario, the angular extent of the flux at 
high energies is smaller than at lower energies. We note 
that decreasing m^ for a fixed m-^ narrows the angular 
extent of the signal, and therefore has little impact on our 
results. We confirmed that for m^ as large as 10 GeV, 
the cross-section limits vary negligibly except for a slight 
weakening of the limit at the lowest end of the m^ range 
considered here. 
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Figs. |6] and [7] show the constraints on ctsi and ctsd as a 
function of m^ , derived from the upper hmits on the solar 
CRE flux obtained in J JInB] For each m^ the CRE flux 
in each energy bin used in this analysis was calculated, 
and the limit on the scattering cross-section was set by 
the energy bin providing the strongest constraint. The 
jagged shape of the curve reflects the transitions between 
the energy bins setting the strongest limit. Models above 
the curves exceed the 95% CL solar CRE flux upper limit 
for the 30° ROI in at least one energy bin. Ref. [l| notes 
that due to the Parker spiral shape of the Sun's magnetic 
field, CREs emitted from the Sun may be deflected in 
such a way as to appear to originate from a source 
displaced by up to 30° from the Sun's position. If we 
instead consider larger ROIs centered on the Sun in order 
to accommodate the expected angular distribution of the 
fiux of a displaced source, the constraints derived on the 
scattering cross-sections would be weakened by ^ 30% 
using the fiux upper limit for the 45° ROI, or by a factor 
of ^ 2 if the 60° ROI fiux upper limit were used. 

The bounds on the scattering cross-sections we derive 
for e^ final states are significantly below the typical 
constraints from direct detection experiments, and so we 
are prompted to examine more closely the validity of our 
assumption of equilibrium. For the limiting values we 
derive on elastic scattering cross-sections, capture and 
annihilation are effectively in equilibrium assuming an 
annihilation cross-section consistent with thermal relic 
dark matter (av) = 3 x 10~^^ cm"^ s~^ for all values of 
the decay length L considered here. In particular, for 
the limiting values of the scattering cross-sections the 
fiux suppression relative to the equilibrium fiux for any 
mass we consider is always less than 3% (following the 
standard calculation implemented in Ref. [36|), and thus 
we work under the assumption of equilibrium, noting that 
there remain uncertainties in the capture rate calculation 
at the level of a factor of a few (e.g., 138}). 

Decays to e^ are generally accompanied by final state 
radiation (FSR), so these scenarios can also be con- 
strained by solar gamma-ray observations. Ref. [39j de- 
rived bounds on the rate of decay to e by requiring that 
the predicted FSR does not exceed the solar gamma-ray 
emission measured by Fermi. However, the constraints 
we obtained on the elastic scattering cross-sections from 
the solar CRE flux correspond to constraints on the an- 
nihilation rate roughly 2-4 orders of magnitude stronger 
than those placed by gamma-ray constraints on FSR. 
The strength of the CRE limits relative to those from 
FSR increases for larger m^. The relative strength of 
the constraints derived from the CRE flux limits com- 
pared to those from the gamma-ray measurements can 
be attributed in part to the fact that the FSR flux pro- 
duced by annihilation to e~^e~ is ~ 2-3 orders of mag- 
nitude smaller than the CRE fiux. FSR emission also 
must compete with a known background gamma-ray fiux 
from the Sun |40j . Furthermore, the FSR constraints 
in Ref. |39j were derived using the preliminary Fermi 
measurement of the solar spectrum which extends only 



to 10 GeV, while this analysis spans CRE energies from 
60 GeV to ~ 1 TeV. Since the FSR photon spectrum 
is harder than the measured solar gamma-ray spectrum, 
the strongest constraints are obtained from the highest 
energy bin in that analysis. Due to the fact that the 
FSR spectrum associated with the DM mass range con- 
sidered in this analysis extends substantially higher than 
10 GeV, the existing FSR constraints are significantly less 
competitive than our CRE constraints. A measurement 
of the solar gamma-ray emission at higher energies could 
likely strengthen the FSR constraints to some extent. 



B. Inelastic dark matter 

We now consider the fiux of e from annihilation 
of DM particles captured by the Sun but with orbits 
which take them outside the surface of the Sun. In a 
standard WIMP scenario, DM particles captured by the 
Sun via elastic scattering quickly undergo subsequent 
scatterings which cause them to settle to the core, and 
hence the fraction of captured DM particles outside the 
surface of the Sun at any given time is negligible 138|. 
However, this is not necessarily the case for inelastic 
dark matter (iDM) [1, IH SI]. This class of models 
has garnered interest recently in light of claims that 
iDM could naturally explain such observations as the 
511 keV hue observed by INTEGRAL/SPI and the 
apparently inconsistent results of DAMA/LIBRA and 
CDMS if the DM scattered inelastically and thereby 
transitioned to an excited state with a slightly heavier 
mass. 

For a DM particle x to scatter inelastically off a 
nucleon N via the process x + -^ ^" X* + ^! the DM 



must have energy E > S{1 



^/itin), where S 



•■X- 



Particles captured by the sun by inelastic 
scattering typically lose enough energy after only a few 
interactions to prevent further energy loss by scattering. 
If the elastic scattering cross-section is sufficiently small 
(^n ^ lO^^^cm^, e.g., Ref. [l|), the captured particles 
will be unable to thermalize and settle to the core, and 
instead will remain on relatively large orbits. As a 
result, the density of captured DM particles outside the 
Sun may not be negligible in an iDM scenario, and the 
annihilation of those particles to e^ could thus produce 
an observable fiux of CREs from the direction of the Sun. 
While it is not necessary for DM to annihilate primarily 
to e* in order to explain the direct detection results 
(since direct detection experiments are not sensitive to 
the dominant annihilation channels), leptophilic iDM is 
strongly motivated since it could provide a consistent 
interpretation of multiple data sets [j, |43 - |44| . 

In the following we will assume that the DM particles 
annihilate at rest and thus the energy of the e^ produced 
in annihilation is well-approximated by i^cRE — "t-x • We 
will further assume that the CREs suffer no significant 
energy losses between production at the surface of the 
Sun and arrival at the detector, and so we expect a mono- 
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energetic flux of CREs in tfiis scenario. 

For simplicity, we assume all annihilations occur at the 
surface of the Sun (as in [l[), since the density of DM 
falls off quickly with distance from the Sun. Naturally, 
e produced in annihilations inside the surface of the 
Sun cannot escape the Sun, and thus do not produce a 
detectable flux. 

The isotropic flux of e* particles from the Sun is 



F = 2 



Ta.o 



At:DI 



(30) 



where Fa. out is the annihilation rate of DM particles 
outside the surface of the Sun. The factor of 2 accounts 
for the fact that 2 CREs are emitted per annihilation 
of a pair of DM particles. However, it is also necessary 
to take into account that CREs produced on the surface 
of the Sun opposite to the Earth are extremely unlikely 
to reach the detector, so we assume the flux of CREs 
observable at the detector is a factor of 2 smaller than 
that given by Eq. 1301 

Following Refs. [4^, |4g], we assume that capture and 
annihilation of particles in this scenario is in equilibrium, 
i.e.. Fa — \Cq, where Fa is the total annihilation rate at 
all radii. We emphasize, however, that due to significant 
uncertainties in the density profile of the captured iDM 
particles, the assumption of equilibrium is less robust in 
this case than in the elastic scattering scenario. Ref. ^M, 
concludes that equilibrium will be attained, but notes 
the sizable uncertainties in this calculation. On the 
other hand, for the limiting cross-sections we determine 
for this scenario, the condition for equilibrium given 
in Ref. |46| for inelastic capture requires a minimum 
annihilation cross-section ranging from more than an 
order of magnitude smaller than for a thermal relic for 
small masses and 5 = 110 keV to a factor of ^ 3 larger 
than thermal for larger masses and 5 — lAQ keV. In light 
of the uncertainties in this calculation, we again work 
under the assumption of equilibrium when deriving limits 
on the scattering cross-section. 

Defining /out as the fraction of captured DM particles 
outside the Sun at a given instant, we have 



TA.out — JoutTA — T^/outC^G- 



(31) 



The capture rate of iDM particles by the Sun Cq was 
calculated by Refs. [43,|4y|- Both studies note that there 
are uncertainties in this calculation at the factor of a few 
level. We use the capture rate as a function of DM mass 
m^ and mass splitting 5 as given in Fig. 2 of Ref. [40| , and 
interpolate the results shown in that figure. The capture 
rates were calculated assuming the following parameters: 
the velocity of the Sun in the DM rest frame u© = 
250 km/s, the DM velocity dispersion v = 250 km/s, the 
local DM density pdm = 0.3GeV/cm'^, and the cross- 
section per nucleon in the elastic limit ctq = 10~^°cm^. 
The relation between the total inelastic scattering cross- 
section and the total elastic scattering cross-section is 
given in Eq. 7 of Ref. [4g| . The capture rate scales linearly 
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FIG. 8: Constraints on iDM model parameters for three 
values of the mass splitting S. Models above the curves 
produce a solar CRE flux that exceeds the 95% CL flux upper 
limit for a 30° ROI centered on the Sun in one or more energy 
bins. 



with pdm and ao, while the dependence on vq and v is 
mild over the mass range of interest {m^ '-^ 100 GeV to ~ 
1 TeV). We note, however, that the constraints obtained 
by direct detection experiments may be more sensitive 
to variations in the assumed velocity distribution of the 
DM particles. 

The parameter /out was calculated by Ref. [l| by 
simulating the capture of DM particles by the Sun via 
inelastic scattering. Here we interpolate the values of /out 
as a function of S shown in Fig. 4 of that work, which 
were calculated for m^ = 1 TeV. Those authors note 
that the dependence on m^ is weak for the mass range 
of interest, thus we adopt the values of /out determined 
by [l| for m^ — 1 TeV for all masses considered. We 
caution that the calculation of /out is subject to severe 
uncertainties, and a detailed study beyond the scope of 
this work is needed to more robustly estimate the value 
of this parameter. In particular, we note that /out varies 
by more than an order of magnitude over the range of S 
values considered in this study, and we therefore stress 
that the calculation of /out introduces uncertainties in 
the derived scattering cross-section limits of at least a 
factor of a few. 

We calculate the flux of CREs from annihilation of DM 
in this scenario as a function of m^ and ctq for three values 
of the parameter 6. We then derive constraints on the 
TO^-fJo parameter space by requiring that the predicted 
flux of each DM model does not exceed the 95% CL upper 
limits on solar CRE fluxes for a 30° ROI centered on the 
Sun, again using the results derived in ijHIBI Since the 
region from which the DM-induced flux originates in this 
scenario is the angular extent of the Sun, the 30° ROI is 
more than sufficient to encompass all of the DM signal. 

The predicted flux is mono-energetic, however the 
flnite energy resolution of the LAT will result in the 
observed events being assigned to more than one energy 
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bin. Since this may have a non-neghgible impact on 
the derived scattering cross-section hmits for DM masses 
near the energy bin edges, we convolve the predicted 
signal from each model with the energy resolution of the 
LAT and calculate its flux in each energy bin used in 
the analysis. We approximate the energy dispersion of 
the LAT as a Gaussian with a given by the half-width 
of the 68% containment window (see Fig. 9 of [7|)- For 
the energy range considered here the energy resolution 
ranges from ^ 5% to ~ 14%. The cross-section limit at 
each mass is obtained from the energy bin providing the 
strongest constraint. 

Fig. [8] shows the constraints from the solar CRE flux 
upper limits on iDM models in the m^-ai^ parameter 
space for three values of 5. Models in the regions above 
the curves exceed the 95% CL flux upper limit for the 30° 
ROI in at least one energy bin. The rounded shape of the 
curves is due to accounting for the energy resolution of 
the LAT. These limits exclude the regions of parameter 
space compatible with the results of DAMA/LIBRA and 
CDMS (in addition to several other direct detection 
experiments) as determined by Ref. [11| for 8 = 120 keV, 
for the range of masses accessible to our analysis {m^ > 
70 GeV) , assuming the dominant annihilation channel is 
e^. Models consistent with both DAMA/LIBRA and 
CDMS at 90% CL exist for values of 6 ranging from ^ 
85keV to ~ 135 keV [ll|; for masses from 70 GcV to 250 
GeV the range of allowed scattering cross-sections is from 
fjo ^ 10-''°cm2 to (To ~ 10-39 cm2 [12]. Although the 
uncertainties in the calculation of the DM fluxes in this 
scenario are significant, we emphasize that constraining 
CTo < 10"^° cm^ is sufficient to exclude the cross-sections 
of models consistent with both data sets. The bounds we 
derive exclude the relevant cross-sections by 1-2 orders of 
magnitude, and hence we conclude that the parameter 
space of models preferred by DAMA/LIBRA can be 
confidently ruled out for m^ > 70 GeV for annihilation 
to e* despite the uncertainties in the fiux calculation. 

This analysis constrains DM models in which the 
primary annihilation channel is to e . We emphasize 
that although other annihilation channels can be probed 
by gamma-ray [sl, HI, H^l or neutrino 0, ill, HI] 
measurements, the upper limits on solar CRE fluxes 
provide a uniquely strong constraint on the e final 
state, which is inaccessible to neutrino telescopes since 
no neutrinos are produced for this annihilation channel. 



V. CONCLUSIONS 

We used a sample of about 1.3 x 10^ CRE events 
with energies above 60 GcV detected by the Fermi LAT 
during its first year of data-taking to search for fiux 
excesses or deficits correlated with the Sun's direction. 



Two analysis approaches were implemented, and neither 
yielded evidence of an enhancement in the CRE flux from 
the direction of the Sun. This result agrees with the 
more general one shown in Ref. [16], where no evidence 
of anisotropics was found in CRE arrival directions above 
60 GeV in the Galactic reference frame. 

Wc derived limits on DM models which generate a 
CRE flux from the Sun's direction for the two scenarios 
discussed in Ref. jl] . In the case of annihilation of DM 
through an intermediate state and subsequent decay to 
e , the upper limits on solar CRE fluxes provide signif- 
icantly stronger constraints on the DM scattering cross- 
section than limits previously derived by constraining the 
FSR emission associated with this decay channel using 
solar gamma-ray measurements. For the iDM scenario, 
the solar CRE flux upper limits exclude the range of mod- 
els which can reconcile the data from DAMA/LIBRA 
and CDMS for m^ > 70 GeV, assuming DM annihilates 
predominantly to e^ . Since direct detection experiments 
are not sensitive to the dominant annihilation channels of 
the DM particles, other data, e.g., solar gamma-ray mea- 
surements and neutrino searches, may be able to further 
constrain these models by excluding regions of parameter 
space for alternative annihilation channels. 



Acknowledgments 

The Fermi LAT Collaboration acknowledges generous 
ongoing support from a number of agencies and insti- 
tutes that have supported both the development and the 
operation of the LAT as well as scientific data analysis. 
These include the National Aeronautics and Space Ad- 
ministration and the Department of Energy in the United 
States, the Commissariat a I'Energie Atomique and the 
Centre National de la Recherche Scientifique / Institut 
National de Physique Nucleaire et de Physique des Par- 
ticules in France, the Agenzia Spaziale Italiana and the 
Istituto Nazionale di Fisica Nucleare in Italy, the Min- 
istry of Education, Culture, Sports, Science and Technol- 
ogy (MEXT), High Energy Accelerator Research Organi- 
zation (KEK) and Japan Aerospace Exploration Agency 
(JAXA) in Japan, and the K. A. Wallenberg Foundation, 
the Swedish Research Council and the Swedish National 
Space Board in Sweden. 

Additional support for science analysis during the 
operations phase is gratefully acknowledged from the 
Istituto Nazionale di Astrofisica in Italy and the Centre 
National d'Etudes Spatiales in France. 

The authors thank Joakim Edsjo for his valuable 
contribution during the preparation of this manuscript. 

JSG thanks J. Beacom, B. Dasgupta, S. Horiuchi, 
D. Malyshev, and I. Yavin for helpful discussions. 



[1] P. Schuster, N. Toro, N. Weiner et al, Phys. Rev. D82, 
115012 (2010). arXiv:0910.1839 [hep-ph]]. 



[2] M. Pospelov, A. Ritz and M. B. Voloshin, Phys. Lett. 



18 



[3] 

[4] 
[5] 
[6] 

[7] 

[8] 

[9] 

[10] 

[11] 
[12] 
[13] 
[14] 
[15] 

[16] 



[17] 
[18] 
[19 

[2o; 

[21 
[22] 

[23] 

[24" 
[25' 



Phys. 



Phys. 
[astro- 



B662 (2008), 53. [arXiv: 071 1.4866 [hep-ph]]. 
N. Arkani-Hamed, D. P. Finkbeiner, T. R. Slatyer 
and N. Weiner, Phys. Rev. D79 (2009), 015014. 
'a^v:0810.0713 [hep-ph]]. 

O. Adrian! et al. [PAMELA Collaboration], Nature 458, 
607 (2009) arXiv:0810.4995 [astro-ph]]. 
J. Chang, J. H. Adams, H. S. Ahn et aL, Nature 456, 
362-365 (2008). 

A. A. Abdo et al. [Fermi LAT Collaboration], 
Rev. Lett. 102, 181101 (2009). |arXiv: 0905.0025 [astro- 
ph.HE]]. 

M. Ackermann et al. [Fermi LAT Collaboration], 
Rev. D82, 092004 (2010). tarXi^^r008:3£ 
ph.HE]]. 

R. Bernabei et al. [DAMA Collaboration], Eur. Phys. J. 
C56, 333-355 (2008). arXiv:0804.2741 [astro-ph]]. 
R. Bernabei, P. Belli, F. Cappella et al., Eur. Phys. J. 
C67, 39-49 (2010). arXiv: 1002. 1028 [astro-ph. G A]]. 
Z. Ahmed et al. [CDMS II Collaboration], Science 327, 
1619-1621 (2010). [arXiv :0912.3592 [astro-ph.CO]]. 
Z. Ahmed et al. [CDMS II CoUaboration], 
arXi v:1012 50781 [astro-ph.CO]] . 

S. Chang, G. D. Kribs, D. Tucker-Smith et al, Phys. 
Rev. D79, 043513 (2009). arXiv:0807.2250 [hep-ph]]. 

D. P. Finkbeiner, T. Lin, N. Weiner, Phys. Rev. D80, 
115008 (2009). arXiv:0906.0002 [astro-ph.CO]]. 

C. Savage, G. Gelmini, P. Gondolo et al, JCAP 0904, 

010 (2009). arXiv:0808.3607 [astro-ph]]. 

W. B. Atwood et al. [Fermi LAT Collaboration], Astro- 

phys. J. 697, 1071-1102 (2009). a rXiv:0902.1089 [astro- 

ph.IM]]. 

M. Ackermann et al. [Fermi LAT Collaboration], Phys. 

Rev. D82, 092003 (2010). ,arXiv:1008.5119. [astro- 

ph.HE]]. 

http : //www . ngdc . noaa . gov/IAGA/vmod/ igr f . html | . 
J. P. Roberts, arXiv: 1005.4668 [astro-ph. HE]. 

E. Orlando and A. W. Strong, Astron. Astrophys. 480, 
847 (2008) arXiv:0801.2178 [astro-ph]]. 

R. Schlickeiser and J. Ruppel, New J. Phys. 12, 033044 
(2010) arXiv:0908.2183 [astro-ph.HE]]. 
E. N. Parker, Astrophys. J. 128, 664 (1958). 
G. B. Rybicki, A. P. Lightman, "Radiative Processes in 
Astrophysics", Wiley- VCH (1979) 

P. Duffett-Smith, "Practical astronomy with your calcu- 
lator". Third Edition, Cambridge University Press (1988) 
http://iau-conim4.jpl.nasa.gov/access2ephs.html] 
K. M. Gorski, E. Hivon, A. J. Banday et al, As- 
trophys. J. 622, 759-771 (2005). astro-ph/0409513]. 
|http : //healpix . jpl . nasa . gov | 



[26' 
[27 

[28' 
[29; 

[30 

[31 

[32; 

[33 

[34; 

[35; 

[36 

[37 
[38; 

[39 

[40 

[41 

[42; 

[43; 

[44 
[45; 

[46 

[47 



G. D'Agostini Nucl. Inst. Meth. A364, 306 (1994). 

G. Cowan, "Statistical data analysis", Oxford University 

Press (1998) 

T. -P. Li, Y. -Q. Ma, Astrophys. J. 272, 317-324 (1983). 

M. Pospelov, A. Ritz, Phys. Lett. B671, 391-397 (2009). 

arXiv:0810.1502 [hep-ph]]. 

I. Cholis, G. Dobler, D. P. Finkbeiner et al, Phys. Rev. 
D80, 123518 (2009). arXiv:0811.3641 [astro-ph]]. 
I. Cholis, D. P. Finkbeiner, L. Goodenough et al, JCAP 
0912, 007 (2009). |arXiv:0810.53M [astro-ph]]. 
M. Kuhlen, D. Malyshev, Phys. Rev. D79 (2009) 
123517. arXiv:0904.3378 [hep-ph]]. 

L. Bergstrom, J. Edsjo, G. Zaharijas, Phys. Rev. Lett. 
103, 031103 (2009). arXiv:0905.0333 [astro-ph.HE]]. 
D. Grasso et al [ FERMI- LAT Collaboration ], As- 
tropart. Phys. 32, 140-151 (2009). arXiv:0905.0636l 
[astro-ph.HE]]. 

M. Cire lh, J. M. Cline, Phys. Rev. D82, 023503 (2010). 
[ar2fiv:1005.1779 [hep-ph]]. 

P. Gondolo, J. Edsjo, P. UUio et al, JCAP 0407 (2004), 
008. astro-ph/0406204 . 

F. Dick and J. W. Norbury, Eur. J. Phys. 30 (2009), 403. 
S. Sivertsson, J. Edsjo, Phys. Rev. D81, 063502 (2010). 

arXiv:0910.0017 [astro-ph.HE]]. 

P. Schuster, N. Toro, I. Yavin, Phys. Rev. D81, 016002 
(2010). arXiv:0910.1602 [hep-ph]]. 

Gigiietto, N., & Fermi-LAT Collaboration 2010, MmSAI, 
81, 99 

D. Tucker-Smith and N. Weiner , Phys. Rev. D64 (2001), 
043502. arXiv:hep-ph/0101138 . 

D. P. Finkbeiner and N. Weiner, Phys. Rev. D76 (2007), 
083519. ;arXiv:astro-ph/0702587 . 

B. BateU, M. Pospelov, A. Ritz et al, Phys. Rev. D81 
(2010) 075004. arXiv:0910.1567 [hep-ph]]. 
I. Cholis, N. Weiner, arXiv:0911.4954 [astro-ph.HE]]. 
S. Nussinov, L. -T. Wang, I. Yavin, JCAP 0908 (2009), 
037. arXiv:0905.1333 [hep-ph]]. 

A. Menon, R. Morris, A. Pierce et al, Phys. Rev. D82 
(2010), 015011. ;arXiv:0905.1847 [hep-ph]]. 
R. W. Atkins, W. Benbow, D. Berley et al, Phys. Rev. 
D70, 083516 (2004). 

The factor a{E) is not exactly equal to 1/100 (that is the 
exact ratio between the sizes of the real and simulated 
event samples), because in the randomization process 
only the overall number of simulated events is fixed, 
but the simulated events in individual energy bins can 
change. However, since the size of the simulated data 
sample is large, a{E) ~ 1/100 in each energy bin. 



